Cosmic Microwave Background bispectrum from primordial magnetic fields on large angular scales 
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Primordial magnetic fields lead to non-Gaussian signals in the Cosmic Microwave Background (CMB) even 
at the lowest order, as magnetic stresses, and the temperature anisotropy they induce, depend quadratically on 
the magnetic field. In contrast, CMB non-Gaussianity due to inflationary scalar perturbations arise only as a 
higher order effect. We propose here a novel probe of stochastic primordial magnetic fields that exploits the 
characteristic CMB non-Gaussianity that they induce. In particular, we compute the CMB bispectrum (b;^ ;^ ) 
induced by stochastic primordial fields on large angular scales. We find a typical value of Zi(Zi + 1)^3 (/a + 
^^^hhh ^ 10^^^, for magnetic fields of strength Bo ^ 3 nano Gauss and with a nearly scale invariant 
magnetic spectrum. Current observational limits on the bispectrum allow us to set upper limits on Bo ~ 35 
nano Gauss, which can be improved by including other magnetically induced contributions to the bispectrum. 

PACS numbers: 98.62.En, 98.7().Vc, 98.80.Cq, 98.80.Jk 
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Magnetic fields are ubiquitous in the universe but their ori- 
gin and evolution is still not fully understood. A popular 
paradigm is that the observed magnetic fields result from the 
dynamo amplification of small seed fields. Astrophysical dy- 
namos have been extensively studied [1] but are not with- 
out their difficulties. Another interesting alternative is that 
the observed large-scale magnetic fields have their origin in 
the early Universe, arising perhaps during inflation or some 
other phase transition (|2tl- If a primordial magnetic field with 
a present-day strength of a few nano Gauss is generated in 
the early universe, it would strongly affect the temperature 
and polarization anisotropies of the Cosmic Microwave Back- 
ground Radiation. Such anisotropies arise both due to met- 
ric perturbations induced by magnetic stresses as well as fluid 
motions due to the Lorentz force llH H ■ We study here a 
new probe of primordial magnetic fields, namely, the possible 
non-Gaussian signals they induce in the CMB. 

Previous work on non-Gaussian signals from a cosmolog- 
ical magnetic field has focused on a homogeneous magnetic 
field l'6'l. Such a field can result in correlations breaking spa- 
tial isotropy and also possibly result in north-south anisotropy 
on the CMB sky. We focus here on stochastic primordial fields 
which are statistically isotropic and homogeneous, and con- 
sider the non-Gaussian signals they induce. 

Non-Gaussianity of the temperature anisotropy of the CMB 
has recently aroused considerable interest in the context of in- 
flationary models. In such models, linearized quantum fluctu- 
ations of the inflaton on sub-Hubble scales lead eventually to 
classical long wavelength curvature perturbations. The Gaus- 
sian statistics of the initial quantum fluctuations lead to Gaus- 
sian statistics for the curvature perturbations and the induced 
CMB anisotropies. Non-Gaussianity can come in these mod- 
els only through higher order effects. The magnetic stresses 
on the other hand, depend quadratically on the field. Hence, 
even for the field with Gaussian probability distribution, the 
magnetic stresses have necessarily a non-Gaussian character 



even at the lowest order This will naturally lead to non- 
Gaussianity in the CMB anisotropies induced by primordial 
fields. Thus bounds on CMB non-Gaussianity can strongly 
constrain primordial magnetic fields, or lead to their detec- 
tion. 

The CMB bispectrum or the 3 -point function of the CMB 
anisotropies, if non-zero, provides an important and useful 
characterization of CMB non-Gaussianity. It contains more 
information, for example, than the single point probability dis- 
tribution function (PDF). We compute here the simplest con- 
tribution to the CMB bispectrum, namely that due to the mag- 
netically induced Sachs-Wolfe type effects. 

The magnetic field is assumed to be a Gaussian random 
field. On galactic scales and above, any velocity induced by 
Lorentz forces is generally so small that it does not lead to 
appreciable distortion of the initial field [7]. So, the magnetic 
field simply redshifts away as B(x, t) = bo(x)/a^, where, 
bo is the magnetic field at the present epoch (i.e. at z = or 
a = 1) 

The magnetically induced Sachs-Wolfe type contribution to 
the temperature anisotropies given in fll^l, which arises on 
large-angular scales, can be expressed as. 
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Here Qsi^) = B^{x,t)/{8T:p-y{t)) = b§(x)/(87rpo), where 
p-y{t) and po are respectively the CMB energy densities at 
times t and at the present epoch. (Note that the AT/T given 
above is obtained for on large-angular scales, just as the usual 
Sachs-Wolfe effect.) Ref. |8!] gives an analytic estimate TZ = 
Rj/2Q as the Sachs-Wolfe contribution, where Rj ^ 0.6 is 
the fractional contribution of radiation energy density towards 
the total energy density of the relativistic component. Further, 
the unit vector n gives the direction of observation and D* 
is the (angular diameter) distance to the surface of last scat- 
ter We have assumed an instantaneous recombination which 
is a good approximation for large angular scales. Note that 
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there could also be additional integrated Sachs-Wolfe (ISW) 
contributions to 7?. |8]. 

The temperature fluctuations can be expressed in 
terms of the spherical harmonics to give AT{n)/T = 
Y.im aimYim{n), whcrc 



aim = 47ri 
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Here VLb {k) is the Fourier transform of VI b (x). Since Q, b (x) 
is quadratic in bo(x), we have tlB{k) — ■^^Tp' / '^^^ ^^i^ + 
s)b*{s) I (Sttpo), where now 6i(fc) is the Fourier transform of 
bo(x). 

A measure of non-Gaussianity in the CMB temperature 
anisotropy is its 3-point correlation function and is called the 



bispectrum, \ 
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. In terms of the a^m's it is given by, 

— <^ Q'll7ni^l2'>T'l2^l3'm3 ^ ■ (3) 
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From Eq.[2]we can express I ^ 
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with Cj23 defined as, 



Ci23 =< ^B{k,)VtB{K)^B{K) > 



Cl23 

(4) 
(5) 



on a vector denotes its unit vector. This result has also been 
obtained in Ref . [,lOi1 . 

We calculate the bispectrum in two limits: (i) the 'equilat- 
eral' case for which the three Z/s are equal, and, (ii) the 'local 
isosceles' case for which ^2 = ^3 ^ ^i- In the former case, 
the presence of [kiD*) in Eq. (01, predominantly picks out 
configurations in wave-number space, in the Vl b bispectrum, 
with all the fc's are almost equal. In the latter case, it picks out 
configurations in k values with fc2 ^ ^3 ^ fci • 

The mode coupling integral can be usefully approximated 
in these two limits following methods discussed in earlier 
works |4, ^]. In case (i), when ki — k2 — k^, we split the 
s-integral into the sub-ranges < s < fci and s > fci . In each 
of these sub-ranges we approximate the mode-coupling inte- 
grand in Eq.|2]by assuming that s <C fci and s 3> fci, respec- 
tively. Similarly, in the second case, when fc2 = fc3 3> fci, 
the s-integral is now split into sub-ranges, < s < fci, 
fci < s < fcs and s > fca, and again in each of these sub- 
ranges, we approximate the mode coupling integrand in Eq.|2] 
by assuming that s <C fci, fci ^ s ^ fcs and s ^ fcs, 
respectively. We also restrict ourselves to spectral indices 
—3 < n < —3/2. In fact blue spectra are strongly con- 
strained by a number of observations, particularly the gravi- 
tational wave limits of Ref. 1 1 1]. For numerical estimates we 
will focus on nearly scale invariant spectra, that is n —3. 
For case (i) we then get 



The magnetic field itself is assumed to be non-helical and 
to have a Gaussian probability distribution. Hence it is com- 
pletely specified by its energy power spectrum say M{k). 
This spectrum is defined by the relation < 6i(k)6*(q) >= 
(27r)35(k - q)Py(k)A/(fc), where Py(k) = (% -fc.fcj/fc^) 
is the projection operator ensuring V • bo = 0. This gives 
< hi >= 2/(dfc/fc)A2(fc), where A2(fc) ^ k'^ M [k) / {2tt^) 
is the power per logarithmic interval in fc space residing in the 
stochastic magnetic field. 

As in fs'], we assume a power-law magnetic spectra, 
M{k) — Ak'"' that has a cutoff at fc = fc^, where fc^ is the 
Alfven-wave damping length-scale We fix A by demand- 
ing that the variance of the magnetic field smoothed over a 
'galactic' scale, fcg — (using a sharp fc-space fil- 

ter) is Bq. This gives, (for n > — 3 and for fc < fcc) 
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The 3-point correlation function of f2s(fe) involves a 6- 
point correlation function of the fields. A long but tedious 
calculation gives C^^g — S{ki + k-z + k^) V'123' where 



123 



TT^ (n + 3)2(7 -n) /fci 
fc"S 2(|n+l|) Ug 



2(|n+l|) 
while for case (ii) we have 
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Here we have defined Va, the Alfven velocity in the radiation 
era as 



Va 



Bo 
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with S_9 = (So/lO-^Gauss). 

We now express the delta function in its integral form 
6{k) = (l/(27r)'^) / cfixexp{ik ■ x), use the spherical wave 
expansion of the exponential terms, substitute it into Eq. (01, 
and integrate over the angular parts of (fci,fe2,fc3,a;). This 
algebra is very similar to what is done for calculating the 
primordial bispectrum II12II . After this algebra we can write 
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f3 = (s ■ s — fca) and 7 = (fci + s ■ s — k^), where the hat 
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and we have introduced the Gaunt integral 
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For case (i), we substitute Eq. (|8]l into Eq. ( fT2] ) for the re- 
duced bispectrum. Similarly for case (ii) we substitute Eq. (|9]l 
into Eq. ( fT2] l. The integrals over ^2 can be immediately done 
using J k^dk2ji.,{k2x)ji^{k2D*) = (7r/2a:2)5(a; - i:**), and 
the delta function makes the x-integral trivial. We are then left 
with integrals over ki and fcs given by 
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Here for case (i) 

C{n) = 
where as for case (ii) 

C{n) 
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Note that for any n we can evaluate the integral in Eq. ( fT4l i 
analytically in terms of Gamma functions. For power law 
spectra, the form of the integrals is the same as the usual 
Sachs-Wolfe term. We focus on the nearly scale invariant 
case, n « —3. Such magnetic spectra are expected to arise 
in inflationary models for primordial magnetic field genera- 
tion that we assume here |2]. Further, let us first consider 
the purely Sachs-Wolfe contribution as in ([§], which gives 
TZ = Rj/20 ^ 0.03. For the equilateral case we then have 



(17) 
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h{h + lMl3 + l)bi^i,i, ~ 2.3x10-23 ( _J_ 



while for the local-isosceles case we get 



h{h+l)k{h+l)bij,i, ~ 1.5x10-^ 

^ " (18) 

Including an ISW contribution can lead to larger signals. We 
remark that for n ^ — 3, the numerical values of li{li + 
1)^3(^3 + l)bi_^i^i^ is higher for the local isosceles case by 
a factor 6.4, compared to the the equilateral case. The 



relative strengths of the bispectrum for the equilateral and 
local case has its origin in the ratio of 1P123 for these two 
cases. It can indeed be seen from equations [8] and |9] that 

'0123(ioca/) = '0123(e9m:i) (96 I " + 1 |)/(| 2n + 3 | (7 - fl)) 

which for n ^ —3 comes out to be ~ 6.4'0i23(eqnii)- 

These values for the reduced bispectrum, should be com- 
pared with a value hih + l)l3{h + ^)bi,i.,i^ 4 x lO-^^WL 
at large angular scales, which arises due to nonlinear terms in 
the gravitational potential, characterized by /nl (cf. |13]). 
Thus the magnetically induced non-Gaussian jignal, due to 
the purely Sachs-Wolfe effect (with TZ as in |S|), is a factor 
of about a few times 10* smaller than the standard signal pre- 
dicted in inflationary models with f^L ^ 1- Conversely if 
observationsconstrain^i(^i + l)?3(/3 + l)&i^Z2i^ < 4x10"^^, 
assuming fp^L < 100 (for example say from WMAP experi- 
ment mi UJ), then we have a limit of Bq < 35 nano Gauss 
on the strength of any primordial field with a nearly scale in- 
variant spectrum. Stronger limits would be obtained if we 
were to use potentially larger contributions from the ISW ef- 
fect. 

We note that the k dependence for the standard primor- 
dial contribution to the bispectrum on large angular scales (for 
scale invariant potential perturbations) is the same as that we 
get here, in both cases (i) and (ii), for magnetically induced 
Sachs- Wolfe effect (for scale-invariant magnetic spectra) (cf. 
lll2[ll3[l ). That is in all these cases + 1)^3 (Z3 + 1)6 ? 
is independent of I's (see Eq. (17) and (18) and Ref. 11311 ). 
Thus the bounds on /jv l got from searching the WMAP data 
for the standard 'local' non-Gaussianity are indeed useful to 
set constraints on Bq. The WMAP limits use the much larger 
range of / values than the range for which the SW contribu- 
tion is important. However the limit on Bq depends only very 
weakly (Bq cx f^J^) on the exact observational limit on Jnl. 
Thus our limits on Bq are expected to be reasonably robust. 

A systematic, statistical study of limits on Bq from the 
CMB power spectrum alone, taking account of all the avail- 
able data and comparing with the sum of scalar, vector and 
tensor modes of the magnetically induced signals, although 
challenging, would be of great interest. Present works deal 
with only a subset of the data or the signals: In typical models 
a field Bq 3 nG for scale invariant spectra, are not excluded 
by any current CMB data (cf. S, H] and the reviews 

UM), while Bq ~ 10 nG would seem to be excluded. In this 
paper we have attempted to probe stochastic primordial mag- 
netic fields, arising from SW effect. The present constraint 
on Bq derived above from the non-Gaussianity limits arises 
from the magnetically induced SW contribution to the bispec- 
trum. Although the SW contribution above is weaker than 
that expected from inflationary contribution, it is of interest as 
this is a novel probe of stochastic primordial magnetic fields. 
This is of value for setting the stage for calculating the full 
contribution by including vector as well tensor modes. Vec- 
tor modes are expected to contribute on small angular scales 
(large I). Tensor modes could contribute over the same range 
of ^-values considered here. Adding these contributions is ex- 
pected to lead to stronger limits on Bq from limits on CMB 
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non-Gaussianity. It is worth noting that if only the magneti- 
cally induced SW contribution is taken into account even the 
limit on Bq from the CMB power spectrum would be much 
weaker. We have focused on SW contribution of scalar modes 
here so as to calculate, as a first cut, the simplest contribu- 
tion to non-Gaussian signals induced by stochastic primordial 
magnetic fields. 

There is a distinct advantage of using bispectrum (as a 
probe for stochastic primordial magnetic fields) over the 
power spectrum. The magnetically induced signal being fun- 
damentally non Gaussian, could be more easily distinguish- 
able in the bispectrum. This is due to the fact that bispectrum 
arising due to magnetic contribution can in principle domi- 
nate that arising from, models of inflation with a small enough 
/nl- On the other hand, the problem with the power spectrum 
arising due to magnetic contribution is that it is generally sub- 
dominant to those arising from inflation generated curvature 
perturbations, for nano Gauss fields and scale invariant spec- 
tra. These will hence require careful analysis to isolate the 
magnetic contribution. The probe using bispectrum developed 
in this paper is free of such limitations. 

In conclusion we have studied a new probe of primor- 
dial magnetic fields, by computing the CMB bispectrum they 
induce on large angular scales. The magnetically induced 
signals on the CMB anisotropics have a necessarily a non- 
Gaussian character, as the temperature fluctuations they in- 
duce depend quadratically on the field strength. The CMB bis- 
pectrum is an important and useful probe of this non-Gaussian 
nature. Earlier work has emphasized the role of the CMB an- 
gular power spectrum in constraining and detecting primordial 
fields. As we have argued here, the bispectrum will provide a 
new and independent handle for constraining such fields. Fur- 
ther numerical work, which also considers all the scalar, vec- 
tor and tensor modes, will help to strengthen the constraints 
derived here. 

Note added: After we had submitted this paper a preprint 
OlTIl appeared which confirms our main conclusions. 
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